Theoretical investigation of the four-layered self-doped high-T^^ superconductors: evidence of pair 

tunneling effect 
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Based on a four-layered self-doped t - J type model and the slave-boson mean-field approach, we study 
theoretically the superconductivity in the electron-doped and hole-doped layers. The neighbor layers are cou- 
pled through both the single electron interlayer hopping and pair tunneling effect. The superconducting gap 
magnitude for the electron-doped band is nearly twice of that of the hole-doped one, which contrasts to our 
previous understanding of the electron-hole asymmetry in high-T^ superconductors but consistent with recent 
angle-resolved-photoemission-spectroscopy experiments in four-layered materials Ba2Ca3Cu408F2. Our results 
propose that the pair tunneling effect is important to examine the multi-layered superconducting materials. 
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The high-Tc superconductivity has not been understood in 
spite that it has been studied intensively for more than twenty 
years fP]. Generally speaking the materials include the CuOa 
layered compound and the Cu02 planes act as the conducting 
planes. The parent compound is Mott insulator and the su- 
perconductivity is realized by doping either holes or electrons 
into the CUO2 planes. 

Recently the discovery of the multilayered material 
Ba2Ca3Cu408F2 (F0234) has attracted much attention HH]. 
This compound has four CUO2 planes per unit cell, which 
can be divided into two inequivalent groups: the outer lay- 
ers which have apical F ions and the inner ones. It is ex- 
pected that this material is a Mott insulator from the valence of 
the Cu ion obtained by the canonical chemical formula, while 
it has been verified by experiments that the compound is in- 
deed superconducting with the critical temperature 60 K. The 
angle-resolved-photoemission-spectroscopy (ARPES) experi- 
ments observed at least two Fermi surfaces |4-6]. These two 
Fermi surface sheets are split along the nodal direction, which 
is different from the bilayer splitting effect observed in Bi- 
2212 material, where the two Fermi surfaces degenerate along 
the nodal direction [7]. In fact, one can conclude from the ar- 
eas of the Fermi surfaces that the two Fermi surfaces are cor- 
responding to electron-doped and hole-doped ones with the 
doping around 0.2 |4]. This result is consistent with previous 
study on multilayered systems that the outer layers have more 
holes than inner layers MS- All of the experiments indicate 
that, due to the presence of the apical F ions, this compound 
is self-doped with the outer layers and inner layers are hole- 
doped and electron-doped, respectively. 

One of the remarkable features of high-T^ superconduct- 
ing materials is the electron-hole asymmetry in the phase di- 
agram, i.e., the superconducting region is much naiTower and 
the transition temperature is much lower in the electron-doped 
materials fioll . It can be explained based on the Fermi sur- 
face topology structure, namely, the hot spots (the crossing 
points of the Fermi surface with the magnetic Brillouin zone 
boundary) for the electron-doped band are close to the nodal 
point while those of the hole-doped one close to the antinodal 
point. Thus when we consider the cZ-wave superconductivity 
meditated by the antiferromagnetic spin fluctuation, the pair- 
ing strength for the hole-doped band is stronger 1 11]. On the 



other hand, the Fermi surface of the hole-doped band is close 
to the antinodal point, which gives rise to the Van Hove singu- 
larities near the Fermi energy. This may be another important 
ingredient for the electron-hole asymmetry |12]. Since the 
electron-doped and hole-doped superconducting planes co- 
exist in the F0234 compound, which may provide an effec- 
tive way to study in detail the electron-hole asymmetry. The 
Fermi surfaces and the corresponding superconducting gaps 
for the electron-doped and hole-doped bands were observed 
recently by the ARPES experiments HIH]. The hot spots 
for the electron-doped Fermi surface are close to the nodal 
point. And those of hole-doped one is close to antinodal point. 
These features are similar to those of single-layered materi- 
als. Thus one may expect that the electron-doped band and 
hole-doped band in F0234 should have similar asymmetric 
behaviors as those observed previously in single-layered com- 
pounds. However, quite surprisedly, the gap magnitude along 
the electron-doped Fermi surface is about twice of that along 
the hole-doped one |4]. This result contrasts to the asymmet- 
ric behavior in single-layer compounds, which suggest that 
there exists other mechanism in F0234 to be responsible for 
the electron-hole asymmetry. 

Theoretically, the t - J model, with no double occupancy 
on each site as an additional constraint, is one of the most 
popular models to study the hole-doped high-T,. supercon- 
ducting system. The electron-doped system can be mapped 
to the hole-doped one by taking the electron-hole transfor- 
mation to satisfy no double occupancy constraint Then 
the electron-doped and hole-doped systems can be treated 
in the same manner For multilayered systems, the inter- 
layered coupling should be taken into account. The Fermi 
surface and other physical properties for the bilayer com- 
pound were reproduced qualitatively by taking into account 
a momentum dependent single particle hopping term f^k 
(cos kjt - cos^j,)^ ifTil [Tsll . Similar bilayer model consider- 
ing one hole-doped layer and one electron-doped layer was 
also proposed to describe the four-layered compound with the 
different layers coupled through single electron hopping (idn 
or interlayer Coulomb repulsion flj, [Tsll . In this way the 
electron-doped band and hole-doped band were obtained suc- 
cessfully while the unusual electron-hole asymmetry cannot 
be explained. The Fermi surface splitting along the nodal 
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direction observed in the experiments cannot be reproduced 
either |16]. Thus so far a reasonable explanation for the 
electron-hole asymmetry and a systematic investigation for 
the multilayered superconducting systems are still lack. 

The motivation of the present work is to fill this void and 
examine this issue theoretically. We note that besides the sin- 
gle electron hopping, another important interlayer coupling is 
the interlayer interaction in the particle -particle channel |19] 
or pair tunneling effect 1I20I1 . It was also shown that the pair 
tunneling term is necessary for a reasonable explanation of the 
penetration depth in YBCO compound |21]. In this paper, the 
interlayer coupling is considered through both single electron 
hopping and Josephson-like pair tunneling process. Based on 
the slave-boson mean-field approach, we show that the exper- 
imental results are explained naturally and the pair tunneling 
term is essential for enhancing the superconductivity. 

We start with a Hamiltonian describing a system with four 
layers per unit cell, 

H^YjHf+H,, (1) 

/ 

where / is the layer index with 1-1,4 representing the hole- 
doped layers and 1-2,3 the electron-doped layers, respec- 
tively. Hp^ is expressed by the t - t' - J model, 

//(') = - y f(')c('>-'c« - y f(')'c«ic®-/z.c.+7« y 
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For the hole-doped layers, cj^ are the creation operators of 

electrons. For the electron-doped ones, c®' are considered 
as the creation operators of holes to satisfy no double oc- 
cupancy condition. Then we have f^' '*^ - -P-^^^ - t, and 

f(l,4)' _ _f(2,3)/ _ ^, 

We use the slave-boson approach to the above in-plane 
Hamiltonian, namely, the physical particle operators c*^ 
are expressed by slave bosons Ijf carrying the charge and 
fermions representing the spin, c*^ — b^l^"^ ff\ At the 
mean-field level, the constraint can be written as bf^^b^'^ + 
Yjo- fia^ fi!^ = 1. In the superconducting state, the pairing op- 
erators are defined as A^'.' = f^P f^-'} - f^^- We assume the 

mean-field order parameters, A[^* = (f^f^i ' fii^n^ " -^^'^ 
(+ are for the bond {ij) along x and y direction, respectively) 
and^*^? = Y,Afi!^^fjtl) = xf- The boson operators bf^ con- 
dense in the superconducting state, bf^ = = V^, where S 
is the doping density . 

The Hi term is the inter-layer coupling, including the single 
particle hoping term Hjs and pair tunneling term H/p, 

Hi = His + Hip (3) 
-^t.Y.f^'f^'^h.c. ' (4) 



Hip - -^Z^f^r-^ZW 

{ij) {ij) 

-Aj]^>^^Uh.c. (5) 
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Then the mean-field Hamiltonian in the momentum space 
can be written as, 

k/o- Ik 

- - cc - Cf^ - h.c.)] . so, (6) 

where 4" = -(2(Jf® + j'-'^xfXcoskj, + cosky) - 
cos k, cos ky-fi'-'\ A^" = A^"*^ = A<i'(cosfe^. - cos;t,,), 
^k* = A^-*^ = A<2'(cosA:j:-cosA;;>,), fj_k = 5t_i_(cos k,-cos ky?)^ /4, 

and so = 2N Y^i J'-'Hxf^ + A("2/2) + 4-NA(2A'-^^A'-^^ + A'-^^^). 
The renormalized gap magnitudes for the hole- and electron- 
doped layers are expressed as. A*" - 7*'^A^'^ + 2/lA*^^ and 
^(2) ^ y(2)^(2) + 2^(A(i> + A^)), respectively. 

Diagonalizing the above Hamiltonian [Eq.(6)], one can ob- 
tain four energy bands, expressed by ^ki/<v - + (i/'fj_k -1- 
V -y/ck - i/'t/k)/2, where i//,v - ±1 and we define e±k = - 
e^^')/2, Ck = 4e^j^ + 5f^|^ and dk - 4e+kfj.k- The energy bands 
for V = -i-l(-l) represent the hole-doped (electron-doped) 
bands, respectively. Then the energy bands in the supercon- 
ducting state can be obtained, /ikrtiv = t -^^ki/iv '^k^^' where 
T = +1, and the layer I - 1(2) for v = H-l(-l), respectively. 

The free-energy is written as, 

F^-2tJ] ln[2 cosh(^)] - (ju^" + //<2))a? + e„. (7) 

kripv 

The mean-field parameters Xq\ and the chemical po- 
tentials f/'^ are calculated from the self-consistent equations 
obtained by dF/dxf = 0, dF/dA^'^ = 0, and dF/dfi^'^ = 
-(1 - 6)N. We set t'-''^^ = -P'^^ = 1, f^-^)' = -f^^') = _o.3, 
= 0.3, and the temperature T = 0.0005. 

The self-doped behavior indicates the charge imbalance be- 
tween the hole-doped layers and electron-doped ones. We de- 
fine the imbalance constant W - -(/!*'' + ijP^^)/2. The system 
is half-filled when the imbalance constant is zero and self- 
doped for non-zero imbalance constant. The imbalance con- 
stant as a function of the doping 6 is plotted in Fig. 1(a). As 
seen, for all the parameters we considered, W equals to zero 
for zero doping. And | W \ increases monotonously as the 
doping increases. 

The gap magnitudes as a function of the doping 5 with- 
out the pair tunneling effect are plotted in Fig. 1(b). As 
t± - A - 0, corresponding to the single-layer superconduc- 
tors, the gap magnitude for the hold-doped system is larger 
than that of electron-doped one. This electron-hole asymme- 
try is consistent with previous experimental and theoretical 
results. When the inter-layer single particle hopping term is 
considered, as shown, the gap magnitudes for the both bands 
change slightly because the Fermi surface topology changes 







/.=0,15 


=}=Q 




).=0.^ - ~ ^ 





5 5 

FIG. 1: (Color online) (a) The charge imbalance versus the doping 
S. (b) The gap magnitude as a function of doping when the pair tun- 
neling effect is absent. The solid lines and dashed lines are for the 
hole-doped and electron-doped cases, respectively, (c) The same as 
panel (b) but with both single particle hopping and the pair tunnel- 
ing effect, (d) The gap ratios (see text) with different pair tunneling 
constants. 



slightly. While the gap magnitude for the hole-doped layer 
is always larger than that of the electron-doped one even for 
a very large single-particle hopping constant, which indicates 
that the single particle hopping term cannot explain the abnor- 
mal electron-hole symmetry in the F0234 sample. 

The gap magnitudes taking into account both the pair tun- 
neling term and the single particle hopping term are plotted 
in Fig. 1(c). As the pair tunneling is weak (A - 0.1), the gap 
magnitude of the hole-doped layer is still larger and the asym- 
metry appears for various doping densities. However, as the 
pair tunneling constant A increases, the asymmetry is reversed, 
namely, for the case of A - 0.15, the gap magnitude for the 
electron-doped layer is larger than that of the hole-doped one. 
This is consistent with the experimental observations |4] and 
our results propose that the pair tunneling effect may account 
for the abnormal electron-hole asymmetry in the F0234 sam- 
ples. 

Our results are summarized in Fig. 1(d). We define the ratio 
77 = A^^V'^'^'*- As seen, for the cases where the pair tunneling 
constant is absent or quite weak, t] is less than one. This asym- 
metry is due to the Fermi surface topology which has been dis- 
cussed intensively in previous studies. As the pair tunneling 
constant increases to 0.15, as seen, the gap ratio reaches near 
2.0. Interestingly, this ratio does not increase as A increases 
further (not shown here). We have checked numerically that 
the ratio is around 1 .6 ~ 2.0 even for very large A. Our results 
agree well with the experimental observation, which gives the 
ratio around 2.0 |4]. 

The origin of the described anomalous asymmetry can be 
elaborated at the mean-field level, i.e., from the mean-field 
Hamiltonian [Eq.(6)], the pair tunneling term can be absorbed 
to the pairing potentials and increases the effective potentials. 
For the electron-doped layer there are two neighbor layers 
thus the effective pairing potential is larger than that of the 
hole-doped one which has only one neighbor layer. As a re- 
sult, the gap magnitude of the electron-doped layer is larger 
than that of the hole-doped one as the pair tunneling constant 
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FIG. 2: (Color online) The normal state Fermi surface (a) and band 
structure (b) with fx = 0.4 and A = 0.15. 



increases. Our results propose that the pair tunneling effect 
should play an important role when one investigates the phys- 
ical properties of the multi-band superconductors. 

We now study the normal state Fermi surface. The pa- 
rameters are fixed as tj_ - 0.4, A = 0.15, and the doping 
6 = 0.2 [corresponding to the imbalance constant W - 0.0735 
displayed in Fig. 1(a)] in the following presented results. As 
seen in Fig. 2(a), there are four sheets of Fermi surface, among 
which two being electron-doped type with hot spots close to 
the nodal point and two hole-doped type with hot spots close 
to the antinodal point. The two distinct bands near the nodal 
direction are observed. Each band splits into two when closer 
to the antinodal direction. Our result of the Fermi surface 
near the nodal direction is consistent with the ARPES exper- 
iments ||J,|5l]. The band splitting near the antinodal direction 
is consistent with LDA prediction |5]. While the Fermi sur- 
face revealed by the experiments is not distinct enough and 
large error bars exist |l4| |5t]. Thus the band splitting effect 
near the antinodal direction has not been resolved by experi- 
ments yet and still needs further verification. The band struc- 
ture along the path (0, 0) -> (tt, 0) -> (tt, tt) -» (0, 0) is plotted 
in Fig. 2(b). As seen, there is one electron-doped band and one 
hole-doped band near (0, 0), respectively. We can also see that 
both bands spilt and the flat regions appear near the antinodal 
point (n, 0) indicating the Van Hove singularities in the den- 
sity of states. For the hole-doped bands the flat regions appear 
near the Fermi energy and the splitting energy is small. For the 
electron-doped ones the flat regions are at the higher energy 
with the splitting energy is larger. The spilt bands combine 
into one as the momentum moves away from (n, 0) towards 
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FIG. 3: (Color online) The density of states in the normal state and 
superconducting state with = 0.4 and A = 0.15. The (blue) solid 
and (red) dashed thin lines are the density of states contributed from 
hole-doped and electron-doped bands, respectively. 

the diagonal direction. 

At last let us study the density of states spectra in the nor- 
mal state and superconducting state, respectively. In the nor- 
mal state, as seen in Fig. 3(a), there are several peaks shown 
up at the negative energy, corresponding to the Van Hove sin- 
gularities as we have mentioned above. For the hole doped 
band, because the band splitting is small from the band struc- 
ture shown in Fig. 2(b), thus there is one peak near the Fermi 
energy. The two peaks at lower energies are contributed by the 



electron-doped band. In the superconducting state, as seen in 
Fig. 3(b), we can observe clearly the superconducting coher- 
ent peaks of the hole bands at the lower energy and those of 
the electron bands at higher energy. For the whole density of 
states, the sharp superconducting coherent peaks contributed 
by the hole-doped bands are seen clearly (labeled by the ar- 
rows). The low energy Van Hove peak shown in the normal 
state is covered up by the superconducting gap. While for the 
electron-doped bands, the superconducting peaks are broad 
and renormalized by the Van Hove peaks. As a result, the 
whole density of states includes two superconducting coherent 
peaks caused by the hole-doped band and several broad peaks 
contributed by the electron-doped bands, including the super- 
conducting peaks and the negative energy Van Hove peaks, 
respectively. 

In summary, we have examined the superconductivity of a 
four-layered self-doped system based on the t - J type model 
and the slave-boson mean-field approach. Both single particle 
hopping and the pair tunneling effect have been taken into ac- 
count to describe the inter-layer coupling. The electron-hole 
asymmetry is obtained for the single-layered case and the re- 
sult is consistent with previous experimental and theoretical 
results. For the case of the four-layered model, similar asym- 
metric behavior is obtained when only the single particle hop- 
ing is assumed. While the asymmetry behavior is reversed 
when the pair tunneling effect is added, namely, the gap mag- 
nitude for the electron-doped band is nearly twice of that for 
hole-doped band. This result is consistent with the experimen- 
tal observation and our results propose that the pair tunneling 
effect should be an important issue to study the superconduc- 
tivity of multi-layered materials. 
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